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Abstract. We present a mechanism that may seed compact stellar objects with stable lumps of quark matter, or strangelets, 
through the self-annihilation of gravitationally accreted WIMPs. We show that dark matter particles with masses above a few 
GeV may provide enough energy in the nuclear medium for quark deconfinement and subsequent strangelet formation. If this 
happens this effect may then trigger a partial or full conversion of the star into a strange star. We set a new limit on the WIMP 
mass in the few-GeV range that seems to be consistent with recent indications in dark matter direct detection experiments. 
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ACCRETING DARK MATTER IN NEUTRON STARS 

There is compelling evidence [1] that as much as 90% of the matter in the Universe is under the form of dark matter 
(DM). A set of very popular candidates for DM are WIMPs (weakly interacting massive particles). This mostly follows 
from what is called the WIMP miracle since cosmological data are consistent with a thermal relic cross section 
< OannihV >~ 10^^^ cm^s^\ WIMPs from the galactic halo can be accreted onto compact massive objects such as 
neutron stars (NS) [2] or white dwarfs. The signature of WIMP self-annihilations (assuming WIMPs are Majorana 
particles) is energy release in the range equivalent to twice the WIMP mass, leading to subsequent particle production 
and/or heating of the star. Current predictions of the WIMP mass span the range from ~ 1 GeV /c^ up to 10 TeV /c^ 
however recent results from direct detection experiments DAMA and CoGeNT seem to favour a light DM particle in 
the range 4-12 GeV/c2 [3]. 

Compact objects such as black holes or NS are considered to be gravitational accretors of DM. NS configurations, 
yielding gravitational mass, M, and radius, R, are computed using an equation of state (EoS) that relates pressure 
to energy density and temperature, p = p{e,T). A typical NS has a mass M« 1.4 Mq, radius R Ri 10 km, surface 
temperature Tsmf ~ lOOkeV and central mass density p « l.Ox lO'^ g/cm^. There is a variety of model predictions 
for NS composition [4] but there is as yet no consensus about the fundamental nature of the matter in NS interiors. 

The accretion rate of WIMPs captured by a typical NS is given by [5] 

m;t;(GeV) Pdm,q 

assuming a WIMP-nucleon interaction cross section a^r^ > lO^'*^ and Pdm.o = 0.3 GeV/cm^. Due to competing effects 
of annihilation and evaporation, the number of accreted WIMPS at time t is obtained by solving the differential 
equation = — r^nnih — Fevap where Fannih is the self-annihilation rate calculated as Fannih — (Canniti*^' )jnldV = 
C^N{t)'^, (cTannihv) is the product of thermally-averaged WIMP self-annihilation cross section and velocity, and 
ny^ = the number density of WIMPs inside the NS, is assumed to be constant. The evaporation rate, Fevap, 

decays exponentially with temperature ^ g-GMm;i;/RT^ negligible with respect to the annihilation rate for NS 

with internal T ^ 0. 1 MeV [6]. With this simplification, the population of WIMPs at time t is given by 

N{t) ^ (^T)tanh(f/T) 



where the time-scale is T = 1 / \fWQ^. For f > > T, when the equilibrium between accretion and annihilation has been 
reached, the number of particles accreted is time-independent, = ^x. 



Assuming a regime when velocities and positions of the WIMPs follow a MaxweU-Boltzmann distribution with 
respect to the centre of the NS, the thermaUzation volume in the compact star has a radius rth = { nj^Gplmj^ )^^^ ~ 

^ (t^)^^^( lo'Vcm^ lOOGeV/c^ -^1/2 ^^^^^ i^^g typical NS conditions as central internal temperature Tc = 10^ 

K, Pc=10'^ g/cm^ and = 1 GeV/c^ - lOTeV/c^, rth « 2 x 10^ - 2 x lO^cm. Then the thermalization volume is 

^th = |?rr3^~ (10-2- 10-8) Vns. 

DARK MATTER ANNIHILATION IN THE NS CORE 

The exact energy released due to WIMP self-annihilation is dependent on the nature of WIMPs and on the output 
product channel. We parametrize this lack of knowledge by assuming a fraction of the annihilation energy is given by 
the efficiency factor, /, typically / ~ 0.01 — 1. The rate of energy released in annihilation processes in the star in the 
thermal equihbrium volume can be expressed as [7], 

^annih = CAN^rHxC^ = f^Mj^C^, (2) 

Taking Pdm=Pdm.o and an efficiency rate / = 0.9, we have that in the range niy^ « 1 — 10^ GeV/c^ ^annih = 
2.74 X 10^5 - lO^^GeV/s. The energy released will convert partly into heat, which can stimulate u,d quark bubble 
formation. Thermal fluctuations may arise in metastable hadronic matter via strong interactions [8]. The timing and 
conditions of these transitions depend strongly on the pressure in the central regions of the star given by the EoS 
of hadronic matter. Then ud matter undergoes non-leptonic weak reactions, such u + d u + s, Xo form drops of 
strange uds matter, which has lower energy as a result of the reduction in Fermi energies through the introduction of 
a new flavor. The stability, among other properties of these drops, commonly caUed strangelets, has been extensively 
studied (see e.g. [9]) and depends on the electrical charge, strangeness fraction and size. These exotic forms of matter 
are currently being searched with the CASTOR in the LHC CMS experiment in the forward rapidity and above earth 
attached to the ISS with the AMS spectrometer [10]. The energy needed to form a stable long-hved strangelet of 
baryonic number A can be calculated from its quark constituents, in a first approximation, using either the MIT 
bag model with shell mode filling or the liquid drop model (for details see e.g. [11]). These calculations show that 
stability conditions point towards large negative electrical charge fractions and high strangeness. In the ideal Fermi- 
gas approximation, the binding energy of a strangelet with baryonic number A composed of massive quarks of flavor i 
{i = u,d,s) is 

E^{Hi,mi,B) = +Ni^ii) +BV (3) 

i 

where A^,- is the quark number in the strangelet of baryon number A = ^(A^u +A^d +Ns), j-ii{np^) is the quark chemical 
potential at baryonic number density ma- The thermodynamical potential for the ith-qaaxk type with mass m; is given 
by Q.i{Hi,mi) = Q-tyV + Q.i^sS + £licC being V = A/ha = ^TtR^ tiie volume, 5 = 47tR^ Hie surface and C = 87tR 
the curvature of the spherical strangelet. Expressions for these potentials can be obtained from [11]. The masses of 
the quarks are taken as niu = 2.55 MeV, m^; = 5.04 MeV, m.y = 104 MeV respectively. In order to consider charged 

strangelets a Coulomb correction term can be added as iscoui = |(t^ + ^^'^ ~ ^^i^<^ ^ — 52^' 

i i 

given the /f /i-quark charge and number density, qi and n,. Then, the energy of a strangelet with baryon number A is 

Estimation of the minimum value of A in a long-hved strangelet is model-dependent but typical minimum A values 
are in the range Aiyij^ « 10 — 600 [12]. Strangelets with 'magic' numbers of quarks are found for bag constant values 
145 < Z?'/^ < 170 MeV. If smaller strangelets than the minimum A are created, they will decay rapidly, however the 
long-lived ones have lifetimes of days. Since this time-scale is larger than the conversion time-scale ?conv ~ lOOi [13], 
it is in principle possible that strange stars could be formed if this conversion is triggered. As an estimate the rate of 
formation would be A^sict = Eannih/E^iet- For / = 0-9^ A = 10, ha = riA.o, = 1 GeV/c^, a number A^siet ~ lO^"' s^^ are 
created. Assuming that in the center of the star, self-annihilations are dominant, the energy deposited in the medium, 
per annihilation, will be Ifm^iP'. This allows us to set a limit on the energy scale, at a given central baryonic density 
Ha. 

Ifm^c^ > EtiMnA),mi,B). (4) 

We now discuss the results, a more detailed explanation can be found in [7]. We have found that WIMP mass shows 
a weak dependence on nuclear physics input, parameterized by B. As the central baryon density increases, the lower 



limit on the WIMP mass decreases below ~ 4GeV/c^. We adopt = Sn^.o as the highest value of the central baryon 
density, although recent pulsar masses constrain central densities to somewhat smaller values. 

In Fig. 1 , we show the WIMP mass to produce a long-hved strangelet as a function of A for an efficiency rate / = 0.9 
and different values of the central density = «yi,o (dash-dotted line), Bn^.o (thick soUd line), 5nA,o (thin soUd line) 
for limiting values of the bag constant. For each pair of curves, the lower line is for B'/"* = 145 MeV and the upper line 
is for B^^^ = 170 MeV. For a typical baryonic density of ~ 3nA,o> a lower limit for the WIMP mass mj(>4 GeV/c^ is 
predicted. 
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FIGURE 1. WIMP mass as a function of long-lived strangelet baryon number A for limiting values of 6^/'* = 140, 170 MeV, and 
wa/wa.o = 1,3,5. 

Since the mass of the WIMP sets a scale of the possible triggering of the strangelet formation we find that as 
more massive WIMPs annihilate, the smaller the efficiency rate becomes that is required to trigger a conversion from 
nucleon to strange quark matter. In summary, through the mechanism of self-annihilation of DM candidates in the 
central regions of a typical NS proposed in this work, we have derived a limit on the mass of DM particle candidate 
using the fact that there is a minimum long-lived strangelet mass needed to trigger a conversion from nuclear matter 
into strange quark matter. In the range m;^ > 4 GeV/c^ the energy released in WIMP self-annihilation is sufficient to 
burn nucleon matter into a long-lived strangelet that may trigger full conversion to a strange star. This work presents 
a scenario compatible with current experimental direct and indirect DM searches and may have relevant astrophysical 
impUcations. 
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